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ABSTRACT 

The dominant linear contribution to cosmic microwave background (CMB) 
fluctuations at small angular scales 1') is a second-order contribution known as the 
Vishniac or Ostriker- Vishniac effect. This effect is caused by the scattering of CMB 
photons off free electrons after the universe has been reionized, and is dominated by 
linear perturbations near the Ry = 2 Mpc/(/ir/0.2) scale in the Cold Dark Matter 
cosmogony. As the reionization of the universe requires that nonlinear objects exist 
on some scale, however, one can compare the scale responsible for reionization to Ry 
and ask if a linear treatment is even feasible in different scenarios of reionization. For 
an f^o = 1 cosmology normalized to cluster abundances, only ~ 65% of the linear 
integral is valid if reionization is due to quasars in halos of mass ~ 10 9 Mq, while 
~ 75% of the integral is valid if reionization was caused by stars in halos of ~ 1O 6 M0. 
In A or open cosmologies, both the redshift of reionization and zy are pushed further 
back, but still only ~ 75% to ~ 85% of the linear integral is valid, independent of the 
ionization scenario. We point out that all odd higher-order moments from Vishniac 
fluctuations are zero while even moments are non-zero, regardless of the gaussianity 
of the density perturbations. This provides a defining characteristic of the Vishniac 
effect that differentiates it from other secondary perturbations and may be helpful in 
separating them. 



Subject headings: cosmic background radiation - cosmology: theory 
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1. Introduction 

While recombination at z ~ 1100 marked the end of ionized hydrogen from the viewpoint of 
a linearly evolving universe, the nonlinear evolution of small-scale perturbations resulted in the 
reionization of the intergalactic medium at much lower redshifts. The fact that quasar spectra 
show an absence of an absorption trough from Lya resonant scattering by neutral H atoms 
distributed diffusely along the line of sight, the Gunn- Peterson effect (Gunn k Peterson 1965), 
means that this reionization must have occurred with a high degree of efficiency before a redshift 
of 5. 

One of the necessary consequences of this reionization is the presence of secondary anisotropics 
in the cosmic microwave background (CMB) due to the scattering of photons off ionized electrons. 
These secondary fluctuations can be divided into two classes: anisotropics due to nonlinear 
structures and linear anisotropies. 

Nonlinear secondary anisotropies are of several types. Some of the more studied of these 
include the scattering of photons off the hot intracluster medium of galaxy clusters (Sunyaev k 
Zel'dovich 1970, 1972; or for more recent treatments see, e.g., Evrard k Henry 1991; Colfrancesco 
et al. 1994; Aghanim et al. 1997), gravitational lensing (see, e.g., Linder 1997; Metcalf k Silk 
1997), the impact of inhomogeneous reionization (Aghanim et al. 1995; Peebles k Juszkiewicz 
1998; Knox, Scoccimarro, k Dodelson 1998), and the Rees-Sciama effect due to the bulk motions 
of collapsing nonlinear structures (see, e.g., Rees k Sciama 1968; Kaiser 1982; Seljak 1996). 

Small-scale linear anisotropies come in fewer flavors. Detailed analyses of linear perturbations 
have uncovered a single dominant effect known as the Vishniac or Ostriker-Vishniac effect (Hu, 
Scott, k Silk 1994; Dodelson k Jubas 1995; Hu k White 1995; Hu k Sugiyama 1996). The level 
of these perturbations has been calculated by several authors (Ostriker k Vishniac 1985; Vishniac 
1987; Jaffe k Kamionkowski 1998, hereafter JK). 

These investigations raise the question of whether a detectable Vishniac effect even exists 
since nonlinear structures must exist on some length scale at the time of secondary scattering 
of CMB photons, as it is only by the formation of nonlinear objects that the universe is able to 
reionize itself. If these scales are comparable to those making the dominant contribution to the 
Vishniac effect, then a linear analysis is inappropriate and a calculation of secondary anisotropies 
must incorporate nonlinear effects. 

In this work we determine the minimum length scale, Ry, which must remain linear in order 
for a linear approach to scattering by ionized regions with varying bulk motions to be accurate 
for the range of angular scales over which one can hope to measure secondary fluctuations. In 
hierarchical scenarios of structure formation, such as the Cold Dark Matter (CDM) model, smaller 
structures assemble at early times, later merging to form larger objects. This allows us to place 
limits on the time between the formation of structures large enough to reionize the universe and 
the time at which Ry becomes nonlinear. At that point, while peculiar velocities of ionized gas 
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continue to be imprinted on the microwave background, the nature of this signature is qualitatively 
different and is best interpreted from another perspective. 

The structure of this work is as follows. In Sec. 2 we describe the Vishniac effect and 
determine the physical length scale on which it depends. In Sec. 3 we compare this to the scale of 
reionizing objects in different reionization scenarios and discuss the applicability of linear theory. 
In Sec. 4 we examine how the Vishniac effect is distinguished from other effects. Conclusions are 
given in Sec. 5, and the various cosmological expressions used throughout are summarized in the 
appendix. 



2. Analysis 

2.1. Approximations 

The Vishniac effect is caused by the scattering of CMB photons by ionized regions with 
varying bulk motions. The temperature fluctuations induced along a line of sight are given by 

AT - /"*() _ _ 

— (Q) = - dta T e- T ^n e (9,t)e-v(9,t), (1) 
1 Jo 

where r(0,i), n e (9,t), and v(0,i) are the optical depth along the line of sight, electron density, 
and bulk velocity, ot is the cross section for Thomson scattering, t is the age of the universe, 
and to is the present age. Following JK, we choose a coordinate system in which 9 represents a 
three-dimensional unit vector along the line of sight, 9 represents a two-dimensional unit vector 
in the plane perpendicular to it, and bold letters represent fully three-dimensional vectors. Thus 
v = (v x ,v y ,v z ), 9 = (01,02,0), and 9 = (9 ll 2 , y/l - 9\ - 9%) « (0i,0 2 ,l), the validly of the 
approximation deriving from the small-scale nature of the effect. Note that n e , v, and r are all 
functions of position, the optical depth being given by r(0,i) = / f *° aTn e (9,t)cdt' , where c is the 
speed of light. 

If we decompose the density field into average and fluctuating components, we obtain, to 
leading order, 

^(0) = _^£ jf 1 dwa l Xe{6 ^ W) (l + S(9w ang , w) - Ar(9w ang ,w))9. V (0^n S ,w)e-^, (2) 

where x e (x, w) is the ionization fraction, tq and At(x, w) are the optical depths due to the average 
and fluctuating density components respectively, a(w) is the scale factor with ao = a(0), £(x, w) 
is the overdensity field defined such that <5(x, w) = p(x,w)/p(w) — 1 where p(x,w) is the density 
field and p(w) the average density as a function of comoving distance, no is the present average 
electron density, and u> a ng is the comoving angular distance, given by Eq. (||). Note that we have 
replaced time by w, the comoving distance defined by dw = cdt/a(t), and rewritten v, and r in 
comoving coordinates, x. Taking the mass fraction of He to be ~ 25%, and approximating helium 
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reionization as simultaneous to that of hydrogen, no = &bPc/ m p x 7/8 = 9.9 x 10~ 6 S7fr/i 2 cm~ 3 , 
where p c is the critical density, m p is the mass of the proton, h is Hubble's constant normalized to 
100 km s -1 Mpc -1 , and is the present baryonic matter density in units of the critical density. 

The Vishniac effect arises by considering the homogeneously-reionizing, low optical depth 
case. In this case x e (x, w) is independent of position and At can be ignored. In practice, 
realistic reionization scenarios lead to low values of optical depth and hence dropping Ar is a 
safe assumption (Hu, Scott, &; Silk 1994). Indeed, the measurement of CMB fluctuations at large 
scales precludes the high degree of damping that would be caused by a high optical depth (Scott 
& White 1994; Hancock et al. 1998). The homogeneity of reionization, however, is a question of 
scales and thus represents a basic assumption on which the Vishniac analysis is based. 

In this limit Eq. @ becomes 

AT f 1 - / ~ - \ 

-jr(G) = ~ J Q dw 9(w)9 • (v^w^ng, w) + q(6»w ang , w)j , (3) 
where q(x, t) = v(x,w)5(x,w) and g(w) is the visibility function 

g(w) = , v, ore T = (1 + z(w)Yx e {w)e , (4) 

ca(w) z c 

with our conventions for the scale factor as in Appendix A. This gives the probability of scattering 
off reionized electrons and is a slowly- varying function of uu. 

Finally, we must approximate both v and q using linear theory. In this case, the density 
contrast at a comoving coordinate x and comoving distance w from the observer is a random 
field with Fourier transform given by (5(k, w) = J (i 3 xexp(— ik. ■ x) <5(x, w). The spatial and time 
dependence of <5(k, w) can be factorized, so 5(k, w) = 5o(k)D(w)/ Dq where <5o(k) = 5(k, 0), D(w) 
is the linear growth factor, given by Eq. (|28|) , and Dq = D(0). The power spectrum is then defined 
by the relation 

<<fo(k)<fo(k')> = (<f (k)<fo*(-k')> = (2^) 3 <5 3 (k + k')P(k), (5) 

where £ 3 (k + k') denotes the three-dimensional Dirac delta function. This completely specifies 
the probability density functional from which 5(k) is drawn in gaussian theories. In the CDM 
cosmogony, P(k) is given by Eq. (|3l|) and Eq. (|33|), and is dependent on the 'shape parameter' 
r, which is given as a function of cosmological parameters by Eq. (|32|) and constrained by 



observations of the galaxy correlation function to be 0.23^q'q34 (Viana and Liddle 1996). The 
overall normalization of P(k) can be fixed by the amplitude of mass fluctuations on the 8 /i" 1 Mpc 



scale as defined in Eq. ( p4|) 



The linear velocity field is simply related to the density field by the continuity equation 
V • v(x) = — a(w)5(x, w) which in Fourier space gives 

_„ N ia(w) D(w) , ~ „ . 

v(k,t) = -A J i-l-^k ( 5o(k), (6) 
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where D(w) denotes the derivative of D with respect to time rather than w, and is given by Eq. 
(|29]). As the velocity always points along the direction of k, only k modes with large values along 
the line of sight can have large peculiar velocities in the 6 direction. But these modes are varying 
with wavelengths much smaller than variations in the window function and cancel when projected 
along the line of sight. It is only the second order (q) term then, that contributes to the integral 
in Eq. (|). 

A simple real-space argument gives a different way of understanding this cancellation. As 
gravitational perturbations to a pressureless fluid can be written in terms of the gradient of a 
scalar potential field (v oc where <j) is the gravitational potential) the curl of the velocity field 
is zero to all orders. Thus a line integral of v along a line of sight is approximately the integral of 
a gradient and is zero except for a small contribution at the end points. From all the terms in Eq. 
only the 5 x § ■ v survives to contribute to -f- . 



2.2. Power Spectrum 

Dropping the velocity term from Eq. ^, we can analytically construct the power spectrum of 
the angular fluctuations in the linear limit. Let us define ^-(k) as the Fourier transform of the 
temperature fluctuations such that ^-(R) = / d 2 9exp(— ik ■ 9) ^-(0), with the angular power 
spectrum denned as -P ang (/-ii)(27r) 2 (5 2 (Ki + K2) = ((-^?-(/?i)^-(/€2))}- At the small angular scales 
appropriate to the Vishniac effect, -P an g(K) is simply related to the usual Cgs used to express CMB 
fluctuations by Ci = -P a ng(^ = ^) (JK). 

Several authors (Vishniac 1987; Kaiser 1992; JK) have derived expressions for the Vishniac 
Cgs. Here we provide a new approach that, unlike other techniques, is easily extended to calculate 
higher-order moments, as is shown in §4. Our method is a simple extension of the usual formalism 
used to calculate single-point moments of the Qf-{Q) distribution. 



r 1 2 

The simplest quantity of this sort is the second moment ( A5 £ (0) ) 



where we use B to 



T 

denote convolution with a beam profile. Given such a profile in Fourier space B(ii), the second 
moment can be calculated as 



( 



d 2 K 

(2vr) 2 ^ iKr/a, " : 



B(KfP ang (K). (7) 



Suppose, however, that instead of asking about the second moment calculated from a single map 
observed with a beam B(k), we instead compute the single-point function as calculated from the 
convolution of two maps, observed by beams B\{k) and B^k). As these profiles are arbitrary, we 
are free to take 



Bi(k) = (2tt) 2 5 2 (k- ki) 



B 2 {H) = (2tt) 2 5 2 (k-k 2 ), 



(8) 
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where S 2 (k) is the two-dimensional delta function. In this case we find 

(^i(O)^(O)) = (2^fP^(K 1 )5\R l + k 2 ), (9) 

recovering the angular power spectrum. Thus the power spectrum of the Vishniac effect can 
be computed from the single-point correlation if we are careful to express the beam profiles in 
sufficient generality. 



Let us consider then 

<^(o)^(o)) = 



T 



T 



dwig(wi) / dw 2 g(w 2 ) 



d 3 k a r d 3 k 



(27r)3y (2tt)3 

i?l(4^1,an g )S 2 (^^ 2 ,ang)e ifca -^ 1+ifc ^ W2 (^(k a , Wl )g z (k b , U ; 2 )), (10) 

where q(k, w) is the Fourier transform of q(x, w). Substituting in the expression for q in terms of 



5: 



Eq. ( |10[) becomes 



ia(w)i>(w)D(w) f A_J o(k%(k _ k ^ 



n 2 



(2tt) 2 



k' 2 



(11) 



(^(0)^(0)) 



dwiG(wi) / dw 2 G(w 2 ) Yl 
Jo i=i 



(2vr) 2 



Bi((fci + k 2 )wi !3 



B 2 ((k 3 + k 4 )w 2)ling )e 



h 2 h 2 



where 



G(w) 



g{w)a(w)D{w)D{w) 



D 2 



(Il^(ki)), (12) 

3 1=1 



(13) 



As G(w) is slowly varying, we can follow Kaiser (1992) in dividing the integrals over comoving 
distance into N statistically independent intervals of width Aw, over each of which G(w) is well 
approximated by a constant. In this case 

N 4 



(^(0)^(0)) = - £ G(w n ) 2 Aw 2 J] 

n=l j=l 



T 



d 3 kj 

(2vr) 3 



)5 2 ((fc 3 + fc 4 ) 



JO 



(k ltZ + k 2tZ )Aw\ . f(k 3jZ + k 4)Z )Aw\ ki <z k 3 



■ 



Jo 



V2. 1,2 \1L 

K l K 3 l=i 



<IFo(ki)>. 



(14) 



As the density fluctuations are taken to be gaussian, we can expand the expectation value of the 
product of overdensities by Wick's theorem, keeping only the terms in which k\ is paired with k% 
or &4. If we then define k' 2 = k\ + k 2 , we find 

,ATBl -(0)^(0)) = Y,G(w n ?Aw 2 j J ^B x {k' 2Wn ^)B 2 {-k' 2 w m ) 



T 



71=1 



1,4 ' 1.2 1 li./ Zv7M2 
1 ill 2 -"-1 1 1 



iWdlki-kxiDjg 



(15) 
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The Bessel function has a width 5k z ~ I /Aw so /s2„z << K/w ang wherever j 2 is appreciable. We 

u 2 ' f 

J ang 



can neglect terms that are smaller by a factor of w 2 k 2 / Aw 2 to obtain 



(^(0)^(0)) 



T 



T 



±G{w n fAw j^j ^B 1 (k' 2 w axis )B 2 (-k' 2 w. 



n=l 

k 2 
K l,z 



k 2 



l.z 



kf fc 2 ||(4o)-ki|| 2 



(2vr) 3 7 (27 

P(h)P(\\(k' 2 ,0) 



angV 



Taking 



Bi(/e) 



Bo(k) = 2-Ka 2 e~ 



(16) 



(17) 



gives the second moment as observed by a beam of gaussian width a, while choosing beam profiles 
as given in Eq. (||) yields the angular power spectrum. In this case 



where 



Py(k) 



— -Pang (ft 







1 dw^^P Y (£/w ang ), 



w 



(18) 



aug 



(2tt)* 



P(fcOP(||((*,0,0) -ki||) 



A:' 2 



1,2 



k' 2 



l.z 



k[ 4 k' 2 \\(k, 0,0) -ki| 



(19) 



which, choosing a coordinate system in which the 2/ axis points along the direction of {k, 0, 0) 
becomes 



U rco pi I 

P v (k) = —J dx J dfiP(kx)P(kyJl - 2xn + x 2 )(l - \x 2 



1 



1 — 2x/jl + x 2 

This is equivalent to the usual expression for the Vishniac power spectrum 



(20) 



Pv(k) 



k 
8^ 



dx I ^ d„P(kx)P(kfi - 2x» + x 2 ) ^ - gjf 



(21) 



as can be seen by rewriting both integrals in rectangular coordinates and applying an origin 
shift. Thus the Cgs are dependent on an integral along the line of sight of a term, Py(k), that is 
independent of redshift and arises from the convolution of the q fields. We note in passing that 
our results are in agreement with Dodelson and Jubas (1995) and are twice the values found in 
JK. 



2.3. Physical Scales 

Having outlined the approximations which are used to calculate this effect and constructed 
the resulting power spectrum of fluctuations, we now examine which physical scales contribute 
most to Py(k). Typically, Eq. (18) is used to calculate the Vishniac effect by integrating over a 
particular matter power spectrum. To study the dependence of the effect on physical scale, we 
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replace the power spectrum that appears in Eq. (21) with P(k)W 2 (kR) where P(k) is the CDM 
power spectrum and W(x) is the spherical top-hat window function, given by Eq. p5|). Here 
we consider a ACDM model in which the current nonrelativistic matter, vacuum, and baryonic 
densities in units of the critical density are Qq = 0.35, = 0.65, and = 0.06 respectively, 
and the "tilt" in the power spectrum as parameterized in Eq. (^Tj) is taken to be flat, n = 1.0. 
If the other parameters are taken to be h =0.65, as = 1.05, T = 0.2 and reionization occurs 
instantaneously and completely at z re = 18, this results in Py(k) and Cg as plotted in Fig. |]. Note 
that these values correspond to ^?s a full order of magnitude smaller than large-scale primary 
anisotropies, pointing out the experimental challenges that must be overcome before secondary 
anisotropies can be measured (see, e.g., Subrahmanyan et al. 1993; Church et al. 1997). For 
reference we also plot the COBE-normalized primary fluctuations as computed by the CMBFAST 
code V2.4.1 (Seljak & Zaldarriaga 1996; Hu et al. 1998; Zaldarriaga & Seljak 1998), for the same 
cosmological model. Note that the position of this line is sensitive only to the Silk damping scale 
and thus is largely independent of the shape of the primordial power spectrum. 

In this figure we see that even when filtered at the 1 Mpc//i scale, the power spectrum is 
appreciably changed, especially at higher wavenumbers. While low I fluctuations are not affected, 
a comparison with the solid line shows that these fluctuations are lost in the CMB primary signal 
and are very difficult to measure. Note that the damping shown in this graph represents not only 
a loss of power in linear theory, but also an increase of nonlinear power. Thus, at the point that 
the 1 Mpc//i scale has become nonlinear, the Vishniac effect is competing with nonlinear effects 
over the range of angular scales in which it is able to be detected. 

By the time the 2 Mpc/h scale becomes nonlinear, however, the peak wavelength of the 
Vishniac effect has been shifted by a factor of 1/2 and ~ 50% of the power of Py(k) has been 
lost. At this point, linear calculations are unlikely to be reliable at measurable I values ^ 4000, 
and a more careful theoretical approach becomes necessary. These length scales are inversely 
proportional to the 'shape parameter' T which is 0.2 for this model. Thus we can conservatively 
fix Rv = 2 Mpc//iFo.2 where To. 2 = T/0.2, as the maximum nonlinear length scale that still allows 
a linear analysis to be appropriate. 



3. Redshift of applicability 

Having determined Ry, we now consider what scenarios of reionization are compatible with 
a Vishniac effect. These scenarios can be roughly divided into two classes: those in which the 
dominant source of ionizing photons is due to stars formed in dwarf galaxies with halo masses 
^ 10 6 M© (Couchman & Rees 1986; Fukugita & Kawakasi 1994; Shapiro, Giroux, & Babul 1994; 
Haiman & Loeb 1997) and models in which reionization occurs due to active nuclei in galaxies 
with halo masses £ 1O 9 M (Efstathiou & Rees 1988; Haehnelt & Rees 1993; Aghanim et al. 1995; 
Haiman & Loeb 1998; Valageas & Silk 1999). See also, however, the issues raised in Madau, 
Haardt, & Rees (1998) and Miralda-Escude, Haehnelt, & Rees (1998), and the more exotic 



-9- 



0.5 
x 0.4 



> 



^ 0.2 
* 0.1 




II IL1± 


1 1 1 1 III 


1 1 1 1 1 III 1 


i i mi i 


- 

nlinf 






r \ \"'\ 

hJ ^i •• 






3 1- 




. i im-H — l 


- 

TTTllll 



0.01 0.1 



1 



k p /h (Mpc) 



10 100 



10 



1 1 



C\2 10" 12 I 



10 



14 



= 1 


Mill V 


MINI 1 1 


1 11111 = 


1 


III 


\\l 1 ll\ \l 1 


I i iKil - 




1000 


10 4 


10 






1 





Fig. 1. — Upper panel: Py(k) from ACDM density fluctuations convolved with top-hat window 
functions of various scales. The solid line is the exact analytic expression, the dotted line 
corresponds to a top-hat filtering scale of R = 1 Mpc/h, the short-dashed to R = 2 Mpc/h, 
and the long dashed to R = 4 Mpc/h. The overall normalization is arbitrary as it is dependent 
on the choice of the scale factor. Lower panel: £ 2 Ci/2-K from density fluctuations convolved with 
top-hat window functions of various scales. Window functions are as in the upper panel. The 
sharply falling solid line shows the primary anisotropics as calculated by CMBFAST. 
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scenarios of reionization described in Scott, Rees, & Sciama (1991) and Adams, Sarkar, & Sciama 
(1998). 

Models in which smaller objects are the most important predict redshifts of reionization 
of z rc ~ 20 while models in which reionization is due to objects on scales ^ 1O 9 M predict 
more modest values of z Te . All models of reionization, however, are constrained by the lack of 
a Gunn-Peterson absorption trough in the spectra of high-redshift quasars, implying that the 
intergalactic medium was highly deficient in neutral hydrogen at redshifts <^ 5. Thus, high-mass 
reionization scenarios can only be successful in cosmologies in which the parameters are such that 
relatively large objects became nonlinear at high redshifts, precisely those models where the linear 
approximation is on the most shaky ground. 

Using our value for Ry from Sec. 2, we are able to determine a redshift of applicability, zy, 
below which the Vishniac approximation is invalid over the measurable range of i scales. The 
number density of objects above a critical over-density, 5 C , is given by Press-Schechter theory 
(Press &: Schechter 1974) as 

dn(M,z) _ [2p(z) S C D 
dM ~ ~\ a(M) 2 D(z) 

where n(M, z) is the number density of collapsed objects per unit mass at a redshift of z, p(z) 
is the comoving density of the universe at a redshift of z, D(z) is the linear growth factor of 
fluctuations, and a(M) is the level of fluctuations on the mass scale corresponding to a sphere 
containing a mass M, which can be computed from Eq. (|34|). Typically, this formula is used to 
determine the number density of virialized halos. In this case 5 c (z) is a weak function of z for 
open models and A models and a fixed value of 3(12-7r) 2 / 3 /20 ~ 1.69 in the = 1 case (Kitayama 
& Suto 1996). 

As a rough rule of thumb we can assume that reionization takes place when the 2a fluctuations 
at the relevant scale have collapsed. In this case, D(z re ) / Dq = 5 c (z rc )/(2\^2a(M)) which is 
~ 0.6/cr(M) in the flat case. We take the linear approximation to be valid up the point at 
which the la scale fluctuations at the Ry scale have reached an overdensity of 1. This gives 
D(z v )/D = \/{V2a{Ry)) « 0.7/a(R v ). 

In Fig. Q we plot both z rc and zy as functions of mass, as it is the mass scale rather than the 
length scale that is most easily identified with different reionization scenarios. We consider three 
cosmologies, representative of parameters that favor both low mass-scale and high mass-scale 
scenarios of reionization. In order to compare with a scenario that is representative of stellar 
reionization, we consider a flat model normalized at the 8 Mpc//i scale (Viana & Liddle 1996). 
Here O = 1.0, A = 0.0, Q b = 0.07, h = 0.5, a 8 = 0.60, and n = 1.0. In this scenario T = 0.44, 
shifting the CDM line and decreasing Ry by a factor of ~ 2. Note that this value of T is 
incompatible with the observed galaxy correlation function. Typical of high-mass reionization 
scenarios, we consider the "concordance model" of Ostriker & Steinhart (1995), which was used by 
Haiman & Loeb (1998a) in their modeling of reionization by quasars. In this case the parameters 



( ^ (M ) (22) 
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are taken to be n = 0.35, A = 0.65, n b = 0.04, h = 0.65, a s = 0.87, T = 0.20, and n = 0.96. 
Finally, we examine an open model with f^o = 0.35, again normalized at 8 Mpc//i. In this case f^A 
= 0.0, Q b = 0.04, h = 0.65, o- 8 = 1.02, T = 0.20, and n = 1.0. 

Let us first consider the = 1 model, represented by the lowest pair of lines. In this model 
o"8 is the lowest of all the cosmologies considered and the perturbations evolve the most quickly. 
The combination of these two effects moves the redshift of applicability down to a value of zy ~ 2, 
and thus one might expect to find an appreciable Vishniac effect. The problem, however, is that 
the low normalization and rapid evolution of perturbations also lowers the collapse redshift of the 
objects responsible for reionization. In this scenario, the 2a, 1O 9 M0 peaks collapse at z re ~ 7. As 
reionization must have occurred with a high degree of efficiency by z = 5, this scenario is only 
marginally consistent with quasar absorption-line observations. Thus a flat cosmology is most 
compatible with low-mass reionization scenarios. If the collapse of 1O 6 M0 halos is responsible 
for reionization, then z re ~ 13 for this model, yielding a larger range of redshifts over which the 
Vishniac effect could be imprinted on the microwave background. This redshift is comparable 
to the revised values calculated in the stellar reionization scenario of Haiman and Loeb (1997; 
1998b), although they consider a somewhat higher range of as values. 

To quantify this further, in Figs. ^ and || we replot Fig. ^, replacing the vertical axis with 
Jo dwG{w) 2 /w^ ng Py{£/w ang ), the contribution to Cg due to bulk motions within a redshift 
of z. We take £ = 4000 in Fig. ||] and £ = 12000 in Fig. ||, and normalize the y axis such that 
the integral is equal to 1 at a redshift of z re (10 6 M(7)). The magnitude of the Vishniac Ci is then 
directly proportional to the vertical width of the gap between J™( Zrea ) dwG{w) 2 /w^ ng Py{£/w ang ), 
and Jo dwG{w) 2 /w^ ng Py{£/w ang ), allowing us to judge the linear and nonlinear contributions 
to Eq. ( |I8| ) in arbitrary scenarios of reionization at a glance. 

From this point of view, reionization by 1O 6 M0 objects results in only a marginal improvement 
the accuracy of a linear treatment. While 35% of the £ = 4000 integral is nonlinear if M re = 1O 9 M0, 
the M re = 10 6 Mq case is still 25% inaccurate. These numbers are somewhat lower in the high-^ 
case, in which 18% of the integral is nonlinear in the low mass case, and 25% in the high 
mass. Note however, that our definition of Ry = 2 Mpc(/?F/0.2) was based on the damping of 
perturbations at £ = 4000. From Fig. [j] we see that perturbations at £ = 12000 are largely damped 
when the matter power spectrum is filtered at the R = 1 Mpc(/ir/0.2) scale, and a more fair 
comparison between Figs. [3| and || would be to shift the vertical lines in Fig. |] to masses lower by 
a factor of 8, yielding much the same numbers as in the £ = 4000 case. 

More typical of high-mass reionization scenarios is the ACDM model represented by the 
dotted lines in Figs. |^-||. In this model, erg is slightly higher than in the flat case and the evolution 
of D(z) is slowed. These effects raise the collapse redshift of 1O 9 M0 peaks to z rc ~ 11, easily 
compatible with Gunn-Peterson tests. Note that our crude estimate of the redshift of reionization 
is almost the same as the redshift of ~ 12 calculated by Haiman & Loeb (1998a) for the same set 
of cosmological parameters, using a more sophisticated Press-Schechter based argument for the 
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Fig. 2. — In each pair of horizontal lines the upper correspond to the ionization redshift 
D(z ie )/DQ = S(z TC ) I (2y/2a(M)) : and the lower lines correspond to the redshift at which a linear 
treatment is no longer valid, D(zv)/Dq = l/(\/2cr(M)). The solid lines correspond to the flat 
model, the dotted lines to the "concordance model," and the dashed lines to the open model, with 
parameters as described in the text. The nearly vertical lines show the mass scales corresponding 
to spheres of radius Ry for each of the models. 
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Fig. 4.— Same as Fig. § but with £ = 12000. The onset of nonlinearity at the R = 1 Mpc(/ir/0.2) 
scale can be estimated by shifting the nearly vertical lines to masses lower by a factor of 8. 
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ionizing flux from quasars. In this scenario zy is also pushed back, although this is lessened by the 
shift of r as compared to the flat case. Thus zy ~ 2.5. 

One might imagine that in this much wider margin of redshifts would lead to an 

accurate linear calculation. Fig. || indicates otherwise. In this case almost 20% of the low-mass 
and over 25% of the high-mass integral comes from redshifts at which Ry is nonlinear. The high-£ 
values are slightly lower, with 12% of the 1O 6 M0 and 17% of the 1O 9 M0 integral taking place 
when Ry is nonlinear, but again these numbers become roughly the same as the I = 4000 case for 
a more fair comparison. 

The most extreme case we consider is the cluster-normalized open model, in which erg is the 
highest and D(z) the most slowly evolving. In this case zy ~ 4, the 1O 9 M and the 10 6 M Q 
schemes reionize at z re ~ 22 and z re ~ 33 respectively. This cosmology yields the largest regime of 
redshift space over which a linear analysis is valid and the most accurate results. Here nonlinear 
Ry scale fluctuations contaminate 15% of the low-mass and 20% of the high-mass C4000 integrals. 
Again these numbers are lower at higher I but roughly the same after accounting for the smaller 
filtering scale of R = 1 Mpc(/ir/0.2). 

As a final check of the validity of our analysis, we construct C^mtered defined as the angular 
power spectrum as given by Eq. ( |T8| ) but replacing P(k) with P(k)W 2 (kR n \(z)) where R n \(z) is 
now the nonlinear length scale at each redshift as in Fig. 2 (D(z)/Dq = 0.7 /a(R n \(z))). In Fig. ^ 
we plot the ratio of Catered to Cg calculated from the unfiltered CDM power spectrum. While 
only a few reionization scenarios are represented in this graph, this nevertheless gives us some feel 
of the accuracy of the linear treatment over different scales, and unlike Figs. 2 - 4, is completely 
independent of our definition of Ry. Here we see that in the range of t values at which the effect 
is most likely to be measured (4000 <^ I <^ 120000) this estimate is in good agreement with the 
accuracies given in the previous figures. Note also that a linear treatment becomes increasingly 
inaccurate with £, and thus measurements of fluctuations at angular scales just below this Silk 
damping scale will be most easy to interpret. 

From these results we can safely conclude that even given our present ignorance as to the 
cosmological parameters, no more than ~ 85% of the contribution to Eq. (18) for measurable t 



values can be calculated by linear theory in a CDM cosmogony. This depends only on the shape 
of the CDM power spectrum and the lack of diffuse La absorption in quasar spectra out to z p> 5. 



4. Is There a Detectable Vishniac Effect, Really? 

At this point, one may raise the objection that our argument is a bit semantic, as there will 
still be scattering due to bulk motions even when linear theory breaks down. Indeed, the Kinetic 
Sunyaev-Zel'dovich effect, which is due to the peculiar velocities of clusters, can be viewed as a 
nonlinear counterpart to the Vishniac effect. Is there not, then, a sort of detectable Vishniac effect 
in cosmological scenarios in which Ry is nonlinear during the epoch of reionization, albeit under a 
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Fig. 5. — Ratio of the contribution to Ct from the linear regime to the total calculated contribution 
as a function of I. The solid lines correspond to the flat model, the dotted lines to the "concordance 
model," and the dashed lines to the open model, and in each pair of lines the upper line corresponds 
to reionization by 10 6 M Q objects and the lower line to reionization by 1O 9 M objects. 
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different name? 



The problem with this point of view is that it overlooks some important physical distinctions 
between these effects. The Vishniac effect is due to the presence of a redshift regime in which 
G(w) is slowly varying and a delicate cancellation takes place due to the lack of curl in the peculiar 
velocity field. As it can be calculated precisely, it provides a unique probe of the reionization 
history of the universe that is not available from measurements of nonlinear fluctuations. 

Furthermore, the Vishniac effect displays a distinct signature of higher-order moments that 
allows it to be distinguished from other contributions. In order to understand why this occurs, let 
us consider the bispectrum B(k\, ^{Ik) 2 5 2 (k\ + k 2 + £3) = (^(«i)^?-(k2)^-(k 3 ))- We can 
apply our 5-function beam approach in the linear regime to calculate this as 
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If G(w) is slowly varying, we can again follow Kaiser (1992) in dividing the integrals over comoving 
distance into N statistically independent intervals of width Aw, 
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As there are an odd number of k z terms, there is no pairing of density fields that does not result 
in an odd k z term. As all the k integrals are even, (^(«i) ^-(k 2 ) t^(«3)) = 0- 

This cancellation can be understood from a more general perspective. B(k\,k 2 ) is generated 
by the expectation value of the triple product of the field q. As q is an isotropic vector field, 
{qi(k\)qj(k2)qk(k-3)) can depend on no vectors other than the k vectors themselves and must 
therefore be proportional to at least one of the them (Monin & Yaglom 1971). This means that 
the i = j = I = z component of this product must be proportional to k± tZ , k 2>z , or k% :Z . Thus the 
odd k z term that results in B(k\, k 2 ) = is due to the isotropy of the density fluctuations. By 
a similar argument, all odd moments of the temperature fluctuations must be zero as these also 
depend on the expectation value of the product of an odd number of q's. Note that this is true 
independent of the gaussianity of the probability distribution functional of S(k). 

This cancellation does not apply to the even moments, however, as an even number of q's 
can be arranged in a way that is not proportional to one of the k vectors. Thus the Vishniac 
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effect is unique in that it is nongaussian independent of the gaussianity of the probability 
distribution functional of £(k), but this nongaussianity is expressed only in the even higher-order 
moments. This alternation of zero-and-nonzero higher-order moments provides a unique signal 
that distinguishes the Vishniac effect from other secondary anisotropies, and provides us with 
the opportunity to use nongaussian statistics as a discriminator between these contributions. 
Note however, that the difficulty of measuring secondary anisotropies may make such an analysis 
difficult to apply in practice. 

5. Conclusions 

Due to the tremendous predictive power of linear theory, comparisons between linear 
predictions and large-scale cosmic microwave background measurements promise to constrain 
cosmological parameters to the order of a few percent (Jungman et al. 1996; Bond, Efstathiou, & 
Tegmark 1997; Zaldarriaga, Spergel, & Seljak 1997). The natural extension of this approach is to 
try to measure small-scale secondary anisotropies and match them to linear predictions to study 
the reionization history of universe. The situation in this case is more muddled, however, as a 
number of nonlinear secondary effects also contribute at these scales. 

The dominant secondary linear anisotropy is a second-order contribution known as the 
Vishniac or Ostriker- Vishniac effect. As this effect can be predicted accurately as a function of 
cosmological parameters, several authors have proposed that its measurement will prove to be a 
sensitive probe of the reionization history of the universe. Reionization occurs by the formation of 
nonlinear structures, however, raising the question of whether a regime of redshift space exists in 
which these objects have collapsed but a linear analysis is still appropriate. 

In this work, we have determined the relevant physical scales that give rise to the Vishniac 
effect in a Cold Dark Matter cosmogony, showing that approximations are already compromised 
when 1 Mpc/(/iro.2o) scales have become nonlinear, and break down when 2 Mpc/(fcTo.2o) dark 
matter halos reach overdensities of 1. The width of the redshift regime over which the effect can 
be imprinted on the CMB is dependent on the cosmological parameters and the reionizing mass 
scale. Schemes in which reionization is due to radiation from active galactic nuclei associated with 
dark matter halos of masses ;> 10 9 M Q are limited by the absence of a Gunn-Peterson absorption 
trough. As reionization must have occurred with a high degree of efficiency before a redshift of 5, 
such models are successful only if one assumes a large value of erg, or considers open models with 
slowly-changing linear growth factors. Both these assumptions push back the redshift at which 
Ry becomes nonlinear, limiting the range over which a linear analysis is appropriate. 

Scenarios in which reionization is due to much smaller objects, such as stars formed in dwarf 
galaxies associated with dark matter halos of masses ;> 1O 6 M0, are able to reionize the universe 
at much larger redshifts even in cosmologies in which as is small and D(z) quickly evolving. This 
represents only a marginal gain however, as the high redshift contribution to the Vishniac integral 
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is roughly proportional to comoving distance, and comoving distances are small at high redshifts. 
Thus low-mass scenarios of reionization are more compatible with a linear analysis not so much 
because they reionize earlier as because they allow Ry to become nonlinear more recently without 
violating Gunn-Peterson limits. 

The Vishniac effect arises from physical processes that are distinct from nonlinear secondary 
anisotropies. Its detection indicates the presence of a redshift regime in which a delicate 
cancellation takes place due to the lack of curl in the peculiar velocity field and slow variations in 
G{w). This leaves a unique signature in the higher-order moments of the temperature fluctuations 
that is absent from its nonlinear counterparts. Furthermore, due to the predictive power of linear 
theory, it represents a sensitive probe of the reionization history not available from measurements 
of nonlinear contributions. 

As with measurements of large angular scale anisotropies, small-scale microwave background 
anisotropy measurements have the potential to uncover much about the history of our universe. 
Also as with large-scale measurements, whether this potential will be realized remains to be seen. 
While the Vishniac effect represents a possible probe of the reionization epoch, the analysis will, as 
always, be more involved than first suggested. Ultimately it will only be through the measurement 
and analysis of small-scale microwave background anisotropies that we will be able to know if 
there is a detectable Vishniac effect. 

I wish to thank Nabila Aghanim, Francois Bouchet, Rychard Bouwens, Andrew Jaffe, Douglas 
Scott, and Naoshi Sugiyama for helpful discussions and am particularly indebted to Joseph Silk, 
whose comments and suggestions have been invaluable during the preparation of this work. I 
thank Uros Seljak and Matias Zaldarriaga for the use of CMBFAST and acknowledge partial 
support by the NSF. 

Appendix 

In this appendix, we provide explicit expressions for the cosmological factors used throughout 
this paper. We allow both Q,q and £l\ to be free. In this case the Friedman equations for the 
evolution of the scale factor of the Universe, a(z) are 

^ = h e(z) = H ^n (i + zf + n A + (i - n - n A )(i + z) 2 , (25) 

and 

°- = H$[n A -n (l + z) 3 /2}, (26) 

where Hq = 100 h km sec -1 Mpc -1 is the Hubble constant, and the overdot denotes a derivative 
with respect to time. 
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We choose the scale factor such that ciqHq = 2c. If we are located at the origin, w = 0, 
then an object at redshift z is at a comoving distance, w(z) = \ J Q Z dz'E^ 1 (z') and at a time 
t(z) = -jjj^ dz'(l + z')~ 1 E~ 1 (z'). Note that this is different than conformal time, defined by 
dr] = dt/a, such that crj(z) = w(oo) — w(z). For any flat universe the angular size distance 
i^ang = w, and in an open universe 



sinh(2uVl — &o — ^a) 
Wang = 2V1-J2o-«a "' 

The growth factor as a function of redshift is 



(27) 



while 



D = a o 5fl Q a (1 + zf 
D a a 2 a [E(z)} 2 D(z) ' 1 j 

The evolution of is given by 

ft(-z) = O (l + z) 3 ^ _2 (^) (30) 

where fio = O(0). 

For the power spectrum, we use 

?7T 2 

P(k) = —6 2 H (k/2) n T 2 (k p Mpc/hr), (31) 

where T(q) is the CDM transfer function, k p = k/ao = kHo/2c is the physical wavenumber, and 
the shape parameter V is defined as (Sugiyama 1995) 

r = n (h/o.5)exp(-n b - n b /n ). (32) 



The factor of 8 in the denominator in Eq. (|3l| ) arises because we are using cloHq = 2c. For the 
transfer function, we use the analytic fit given by Bardeen et al. (1986) for the CDM cosmogony, 

rni \ = ln(l + 2.34g)/(2.34g) 

{q> [l + 3.89g+ (16.1g) 2 + (5.46g) 3 + (6.71g) 4 ] 1 /4- 1 > 

The normalization factor 5h can be fixed by specifying er(8Mpc//t), where the variance of the 
mass enclosed in a sphere of radius R is given by 

1 r°° 

<j 2 (R) = — / k 2 dkP{k)W 2 {k p R). (34) 
2tt z Jo 



Here W(x) is the spherical top-hat window function, defined in Fourier space as 



W(x) 



sin(x) cos(x) 



x 3 



(35) 
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